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Abstract. In the 1960s, the technique of using cluster expansion bounds in order to 
achieve bounds on the virial expansion was developed by Lebowitz and Penrose (1964) 
and Ruelle (1969). This technique is generalised to more recent cluster expansion bounds 
by Poghosyan and Ueltschi (2009), which are related to the work of Procacci (2007) and 
the tree-graph identity, detailed by Brydges (1986). The bounds achieved by Lebowitz 
and Penrose can also be sharpened by doing the actual optimisation and achieving ex- 
pressions in terms of the Lambert W-function. The different bound from the cluster 
expansion shows some improvements for bounds on the convergence of the virial expan- 
sion in the case of positive potentials, which are allowed to have a hard core. 



1. Introduction 

In this article, further bounds on the virial coefficients and consequently the radius of 
convergence of the virial expansion are developed in a similar way to how they are achieved 
in the paper by Lebowitz and Penrose [Leb64] . that is via cluster expansion bounds. These 
techniques are shown in general and they rely upon the Lambert W-function. They give 
bounds which improve upon those given by Lebowitz and Penrose and they are applied 
to a more recent version of cluster expansion bounds, detailed by Poghosyan and Ueltschi 
Pog09]. The new bounds display an improvement for small temperatures (large /3) and 
potentials with a hard core and sufficiently attractive tail. There is no difference for purely 
repulsive potentials. 

1.1. Background to the Paper. The bounds on the virial coefficients and consequently, 
the radius of convergence of the virial expansion, were mainly achieved in the 1960s by 
Lebowitz and Penrose lLeb64] , Groeneveld }Gro62] and Ruelle [Ruc69 . These bounds have 
not been improved upon in the literature, even though our understanding of cluster and 
virial expansions has improved since then. Moreover, further, more refined, techniques have 
been developed, which mainly focus on tree-graph identities. Recent work by Pulvirenti 
and Tsagkarogiannis jPu!12j and Morais and Procacci |Morl3j have made some gains in 
understanding bounds on coeffcients, through the use of Canonical Ensemble calculations, 
which naturally use the parameter of density, rather than the fugacity, which is the parameter 
in the Grand Canonical Ensemble. This saves needing to invert the fugacity as a function of 
density or using Lagrange inversion ideas. Furthermore, the recent paper by Jansen |Janl2j 
raised important questions about achieving the low temperature asymptotics, which are of 
the form: e _,9B , rather than e~ 2/3B , which encourages a search for an improvement to the 
current bounds. 
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2. Main Results 

The main results of this paper rely on the particle model of the classical gas. The main 
parameters in the model are: the inverse temperature j3 — t^, where k is Boltzmann's 
constant and T is the temperature; the fugacity z\ the pressure P; and the density p. 

We start with the Grand Canonical Partition function for a classical gas: 

oo 

E(z) = z N* n (2-1) 

N=0 

Where is the TV-particle configuration integral: 

1 Jl/f , \ -f> E 



ZN = mUXL ddxi r (2 ' 2) 

We have the relationship j3P = InS, which gives us an expansion for pressure. The power 
series for pressure in terms of fugacity z, may be written as: 

(3P = Y,b n z n (2.3) 

n>l 

Using the formula p = zJ^((3P), we have the corresponding fugacity expansion for density: 

p=J2 nb n* n (2-4) 

n>l 

We see that the power series expansion for p has a zero constant term and non zero z 
term, where we have that b\ ^ and usually we set it to 1. It is usual to have a non zero 
z coefficient in the cluster expansion. This means that it is possible to invert the p — z 
relationship, in order to obtain an expansion for z(p). This can be substituted into (2.3 1 in 
order to get a power series for P in terms of p: 

p = J2 c ^p n ( 2 - 5 ) 

n>l 

The assumptions for the initial cluster expansion bounds used in this article are: 

Assumption 1 (Potential). The N -particle interaction potential: Un{x\, ■ • ■ ,xn) may be 
written as the sum of pair-potentials: 

U N (xi,- ■ ■ ,x N ) = ^ ®{ x h x i) (2-6) 

Furthermore, we assume that the pair potentials $(xi,Xj) are central, that is they only 
depend on the distance from Xi to Xj . 

Assumption 2 (Stability). The potential energy is assumed to be stable, that is, there is a 
B > 0, such that for every N and (x±, • • • , xn) £ M. Nd , we have: 

U( Xl ,--- ,x N )= ®{xi,Xj)>-BN (2.7) 

l<i<j<N 

where d is the dimension of our system. 

Definition 1 (C(/3) and R{(5) - 'temperedness'). We have two main functions of f3, which 
play an important role in the cluster expansion bounds: 

COS) := / e-W '*) - 1 d d x (2.8) 
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R(13) :-- 



V 



\ 



\r d + p J \$(0,y)\d d y 

\y\>r 



(2.9) 



If the expression (2.8) is finite then the potential $ is called 'tempered', which we assume 



for bounds involving C(j3) 

The r in (2.9) represents the radius of the hard-core interaction. |B| is the surface area 
of a d-dimensional sphere. 

We write lZ\f lr for the radius of convergence of the virial expansion. 

Definition 2 (Lambert W-function). We denote by W(z), the Lambert W-function with 
domain R + := {x G Klimx > 0} and range M + . It is the solution to: 



W{z)e 



W(z) _ 



(2.10) 



Further notes on the Lambert W-function can be found in |Cor97| 
Theorem 2.1 (General Virial Bounds). Assuming cluster coefficient bounds of the form: 



\nb n \ < a- 



(2.11) 



where a and b are non-negative functions of inverse temperature j3, we have the virial coef- 
ficient bounds: 

n-i 




W 



( eab \ 
1+ab J 



^(ff^)- 1 )' 

which gives the lower bound on the radius of convergence as: 



(2.12) 



K 

Vir ^ a 



1 



(2.13) 



If we apply this general theorem, which is derived in Section [4j to two specific bounds we 
have for cluster expansions, we achieve: 



Corollary 2.2 (Improved Lebowitz- Penrose). The cluster expansion bounds: 



\nb n \ < 

give the bound for virial coefficients: 

( 



2PB{n-l) 



n-1 



C n < 



I3- 1 



C(/3)e 



4/3 B _ 



n-1 



V 



l + e 2fll 



(2.14) 



(2.15) 



and the bound for the radius of convergence: 



n Vir >C{[3)- l e 



-1-4,35 



w 



( l+e^s J 



(2.16) 
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This is shown in Section [5] 

For purely hard-core interactions B = and the radius of convergence satisfies: 

Kvir > C^r 1 ^ W ^ (2.17) 

This is precisely the same as what is obtained by Lebowitz-Penrose. 
Corollary 2.3 (Alternative Bounds). For cluster expansion bounds: 

n-l 

\nb n \ < —RiPT^e^ 3 (2.18) 
n! 

we have the bound for virial coefficients as: 

\Cn\ < ^RW)^ 1 - , K \ \ 2 (2.19) 



n \ I W / e' 3B + 1 



1+eP' 



- 1 



TTius giving the radius of convergence as: 

2 



(/ 8B + 1 



1 

K Vir > R(py L - ^ — ' , y (2.20) 

Remark 1. T"/ie bounds for the cluster coefficients used in Corollary \2.S\ are obtained in 
|Rue69| and |Leb64| . The bounds for the cluster coefficients used in Corollary 2.3 are ob- 
tained in |Pog09|. 

Having a new separate bound is only a good idea, if it is an improvement in certain cases. 
In Section [7] we understand how the two different bounds compare and conclude: 



Proposition 2.4 (Comparison of Bounds). The bound in Corollary 2.3 is better than that 
in Cor pilar y \2.2\ precisely when 

1.6R(f3) < C{f3) (2.21) 

This is likely to happen at low temperatures (large /3) and for potentials with a hard core 
and attractive tail. The attractive element of the potential is important, since the improve- 
ments are made corresponding to the stability parameter B. 



3. Obtaining the Relationship 

Lagrange Inversion techniques are the usual approach to obtaining the virial coefficients 
and this has a nice representation as a contour integral. If we want to compute the coefficient 
c„, we can take the contour integral (around 0) of ^ divided by np n . This gives us the 
formula: 
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We can manipulate this equation to get it in terms of the z-variable, since we know about 
the cluster expansion already. 

1 r 

9z -dz 



dz 

(3.2) 



2iri Jqi np' 
fi-l r .dp 1 '. 



2niri J c , zp n 

z 3 " 1 I dz 



2mri Jqi zp n 1 

This (implicit) relationship between the virial and the cluster coefficients is our starting 
point. The idea is that we can bound such an integral quite easily when we consider bounds 
on |/3 1, when we write it in terms of the fugacity z as in (2.4). 

Indeed this method is the one used by Ruelle Rue69j and Lebowitz and Penrose |Leb64j . 

4. General Derivation: Derivation for Theorem 12.11 



We wish to obtain bounds on virial coefficients using those which come from cluster 
coefficients, in order to gain an estimate of the lower bound for the radius of convergence of 
the virial expansion. The derivation loosely follows the method used by Ruelle in |Rue69j 
and that of Lebowitz and Penrose in [Leb64 . The starting point is the cluster expansion 
bound, which is given in the form: 

n n-l 

\nb n \ < a^b n (4.1) 

Where a and b are positive functions of temperature defined by the particular bound we use. 

The starting point is the relationship in (3.2), (implicitly) defining virial coefficients in 
terms of cluster coefficients: 

c "=2^i^T (4 - 2) 

In order to get an upper bound on |c„|, we need to bound \p\ from below. To do this we use 
the bound: 

oo 

\p-z\<Y,\nb n \\z\ n (4.3) 

n=2 

We then substitute in for the upper bounds we have on the cluster coefficients: 

n=2 



We define the function: f(x) := n , x n and write (4.4) conveniently as: 

\p-z\<a(f(b\z\)-b\z\) (4.5) 
Using the reverse triangle inequality, we obtain: 

\p\>\z\(l + ab)-af(b\z\) (4.6) 
Let b\z\ — se~ s and observe that, from the assumed generic bound on the cluster expansion, 

oo 

b\z\ < e~ l in order for ^2 nb n z n to converge. se~ s is an increasing function on (0, 1) and it 
takes values in (0, e _1 ) as required. We then have the bound in terms of s: 

\p\ > b^se-^l + ab) - af(se- s ) (4.7) 
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We have chosen the function se~ s since it is the inverse function of /. This can be understood 
from Lagrange Inversion (shown in appendix) and is a result in |Cor97j . We thus have the 
expression: 

\p\ > 6 _1 se~ a (l + a6) -as (4.8) 

We thus seek to maximise the right hand side to get the best possible bound. We notice 
that in the range s £ (0, 1) we have a zero at s = and another when b~ l e~ s (l + ab) — a = 
i.e. when e s = 1 + 4-, so at s = In (l + X). It is positive for s e (0,ln (l + 4))- We seek 
the value of s to maximise this function in this range. 

Remark 2. This approach to estimating a lower bound for the radius of convergence for 
the virial expansion, takes a value of \z\ some distance away from its maximal value for con- 
vergent cluster expansions. We expect density to increase with fugacity and so the maximal 
density for which the cluster expansions remain convergent appears to be greater than that 
for the virial expansion. This is not true in general and is a weakness of this approach. 

If we define: 

r(s) := s(e _a (l + a&)& _1 -a) (4.9) 
and take the derivative to search for an extremum in this range. 

r'(s) = (e~ s (l + ab)^ 1 - a)- se~ s (l + ab^ 1 (4.10) 



if we find when r'(s) = in (4.101, then we solve: 



(1-*--^ (4.11) 
Substituting 7 = 1 — s, we get the equation for 7: 

The Lambert W- function, as is explained in |Cor97] , is the inverse of 7e 7 and so we can write 
(4.9 1 in terms of 7 and substitute 7 for W(/x), where /i := o^g - If we write r(-/(s)) :— r(s). 
We thus get: 

^(7) = C 1 - l)( e ' y l+ ^ b - a ) 
eb 

l + ab 1 + ab 

= a(e ' , 76 ' 1 + 7) 

eab 

= o(i-2 + 7 ) 

7 

(W{H) -if 



eat 



(4.13) 



Where we use p~12) to cancel: 7e 7 ^ = 1 an d = I. 

We evaluate the integral (4.2) along the contour described by the circle \z\ = constant, 
where the constant is determined by the manipulations above. This leaves us with the 
integral: 

2mr Jc \z\\p\ n 1 
This gives us bounds on the coefficients c„ as: 

m < ^ u-j^r (4.i 5 ) 



n V (W(ji)-1) 
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This gives us the radius of convergence of the virial expansion satisfying: 

(Wfa) - l) 2 

Kvir - a wfr) (416) 

5. Derivation of Corollary [2J2] and comparisons with previous results 



5.1. Derivation of Corollary 2.2 , We have the bound on cluster expansions from |Rue69j 



as: 



\nb n \ < ^e^^CW)- 1 (5.1) 

Where, the parameter C(/3) < oo defines a tempered potential, which is defined in ( |2.8[ ) , and 
B is the stability parameter from Assumption [2] This corresponds to setting the parameters 
in the previous section to: 



a 



This gives: 



= C{0)- x e-^ B (5.2) 

b = e 2f3B C{P) (5.3) 

ab = e~^ B (5.4) 
_ eab _ e 

M ~ 1 + ab ~ l + e^B W 



\c n \<^(c(^ B | (5.6) 

(^( TT iW)-l) 2 



n 



with the radius of convergence satisfying: 



ttvfc>C(fl-V ^^ B ) ^ (5-7) 



( l+e 2 ^ B ) 



5.2. Morais-Procacci Bound. This is an alternative derivation of the Moraid-Procacci 
bounds using the method outlined in Section [4] 

In the recent paper [Morl3 , Morais and Procacci obtain a bound for the virial coefficients 
via the Canonical Ensemble, using polymer expansion methods and the bounds for cluster 
coefficients, in the form: 

\nb n \ < " e^-aJcGfJ)- 1 (5.8) 
n! 

The values a = C(/3)" 1 e~ 4/3B ; b = e 2 ^ B C(i; and thus ab = e" 2/3S are substituted into 



(4.9 1, so that we get: 

r{s) = s(e- s (l + e-^ B )C{P)- l e- 2 ^ B - C^)"^" 4 ^) (5.9) 
If we let u — e 2l3B , then our equation becomes: 

r ( s ) = C(p)- 1 -s(e- s (l + -)--) (5.10) 
u u u 

We reiterate the remark that s € (0, ln(l + u)) so that we change variables (monotonically) 
to s = ln(l + u(l — e~ a )), so that: 
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r(a) 



1 

C(f3)u 
1 



ln(l + it(l-e- Q )) 
ln(l + u(l - er a )) 



u+1 
u(l + u(l - e~ a )) 



C{P)u e a {l + u{l - e~ a )) 
This then gives us the bound for \p\ as: 



\p\ > cwr 1 



ln(l + it(l - e~ a )) 



(5.11) 
(5.12) 

(5.13) 



ae(o,oo) ue a (l + u(l - e -a )) 
If we then follow the same argument before, we get the bound for the virial coefficient as 

?; 



< 



(5.14) 



i -7-/ \ ln(l+u(l — 

where F{u) = max ae(0iOo) MeQ(1+i ; (1 _ 



")) 
"))■ 



The main improvement in the paper comes from better techniques of approximating virial 
coefficients. Using techniques involving Canonical Ensemble calculations for free energy and 
its relationship as the Legendre transform of pressure, they obtain slightly better bounds on 
the virial coefficients, although still under the proviso that \p\ < p*, where p* is the radius 
of convergence in (5.7 1. The improved estimates for the coefficients seem to imply that there 
may be a way of extending the radius of convergence, at least for temperatures for which 
the bound on the coefficient is an improvement. They obtain bounds on the coefficient of 
p k+1 in the free energy as: 



(5.15) 



where a*p is the optimal a from (5.13). This leads to the (better) asymptotic bound for the 
virial coefficients as: 



K 



0.24026 



) 



(5.16) 



where if is a constant. 

5.3. Lebowitz-Penrose. We again use the same bounds on the cluster expansion and 
follow the method outlined in }Leb64j and starting with (5.10), we arrange our expression 
r(s) into the form: 

K*H^^((i + ") 2 e - s -(i 



We use the identity: ^ L 
r{s) = 



1 + u u 2 

(l+u) 2 1+u 



(5.17) 



1 + u 



1 



, to rewrite it as: 
-u 2 ( l + u ) 



We make the change of variables v = and define g(w) 



■s 

l-e" 



and obtain: 



CQ8)- 1 , 
r{s) — (v 



v 2 g 



1 



(5.18) 



(5.19) 



1 + u 

We note that g'(w) - ^"-^-^) = ^+^-1 

If we use the inequality e w > 1 + w, then 1 > (1 + w)e~ w and so > (1 + w)e~ 
hence g is decreasing. We note that u £ (1, 00) and so e (^, 1). Since g is decreasing, we 
get maximum value at <?(|), which would give a minimum value for the expression in brackets 



1 and 
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which gives a bound uniform in u. This is an extra approximation in the derivation of the 
virial coefficient bounds, which is not made in the general derivation. This is equivalent to 
realising that the g- function is dominated by its contribution t u = 1, that is when /? = 0, 
so at high temperatures. 

We then seek to maximise: f(v) — v — v 2 g(\v). Upon differentiation we get: 

f{v) = l-2vg{\v)-\v 2 g'{\v) 

= 1 - 4(1 - e~i v ) - 2(^v + l)e-3" + 2 

= -l + 2e~i v -ve~i v (5.20) 
When we set this to zero we get: 

l = (2-v)e~i v (5.21) 
changing parameters to S = 1 — ^v, we have: 

\ = Se 5 (5.22) 
and so using the Lambert W-function again we get: 

S = W (D ( 5 - 23 ) 

In our original expression: 

f(S) = 2 - 26 - (2 - 2S) 2 g(l - S) 

/I - e 6 - 1 

= 2(l- ( 5)-4(l-5) 2 ( i T ^ 

= (l-^)(2-4(l-e 5 e- 1 )) 

= 2(-l + 6+ 1 -l) 
o 

= 2- — ^ — (5 24) 



This therefore gives us the bound 



2CW)- 1 W(§) - I) 2 
H " l + u W(l) (5 ' 25j 



Leading to the coefficient bound of: 



and the lower bound on the radius of convergence as: 



C (0)-i 1 2 ( ^ ( ^' 1)2 f5 27) 
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5.4. Comparison of the Bounds. If we define: 



(^(l+iw)-l) s 



(5.28) 



jy 2 (^(l)-i) s 



(5.29) 



In order to concentrate on the factor over which the two bounds (5.271 and (5.7) differ, we 



need only concentrate on r% and r-i above. In Figure [Tj we see that the optimised bound 
shows a slight improvement over the Lebowitz-Penrose bound. Furthermore, considering 
the quotient ^ in Figure |2j we see that the optimised bound is 1.25 times better for the 
low temperature limit, whereas at high temperatures the two bounds are approximately the 
same. This is explained by emphasising that the approximation on the g-function is its exact 
value at zero /3 or high temperature. 



Optimised Bound 



Lebowitz-Penrose 



Figure 1. Comparison of the Optimised Bound (ri) (5.28) with the 
Lebowitz-Penrose Bound (ra) (5.29) 
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FIGURE 2. Quotient of the Optimised Bound (n) ( |5.28 1 and Lebowitz Pen- 
rose Bound (r 2 ) (5.29) 



6. The Alternative Bound: Corollary 12.3 



There are alternative cluster coefficient bounds relating to Tree-Graph-Identities found 
in |Bry86| , |Pro07| and [Pog09 1 . The bounds are given by: 



\nb n \ < -—RiPT^e^ 13 (6.1) 
n\ 

where 

R{p)= \n\r d + (3 J |<J>(0,y)|d d J (6.2) 

\ l»l>r / 

where |B| denotes the surface area of the unit sphere in d-dimensions. In this case our 
parameters are: a = b = i?(/3)e' 9S ; ab = e* 35 ; and fi 



The general bound for \p\ is (from (4.13)) 



l+e* 



-PB+1 



Which gives the bound on the coefficients as: 



P\ < R VP) „Z e vB+i : (6-3) 
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and the lower bound on the radius of convergence is: 



rw)- 1 ' Kl T ' x ' (6.5) 

W I e 



7. Comparison of the Separate Bound: Proposition 12.4 



This approach would be better for potentials where (R(/3) < C(P), for some C, representing 
the quotient of the coefficients fy- , where: 



2 

-1' 



^ > (7-1) 



( l+e 2 f B ) 



y J.-1-e- 

2 



1 



/ 2 := V V 7 J +1 , ' (7.2) 

/2 is actually smaller than /i , as detailed in Figure [3] This would only be an improvement 
if the i?(/3) factor compensates. In the case of a hard-core contribution, both R{j3) and 
C{0) are the same, and so are the factors. The use of the integrand j3\U(y)\ rather than 
| e -^f (j/j — 1 1 is better for potentials with attractive parts. For (3 — 0, the coefficients /i and 
/2 are precisely the same. 



The comparison between the optimised bound in (7.1 ) and the bound achieved from these 
different cluster coefficient bounds based on tree graph identities (fa) is shown in Figure [3j 
This gives us that as soon as 1.6R(/3) < C(j3) this other version of bounds is better. This is 
due to the fact that 

C\P) 1 _ /l „v 

Rj/3) h (? - 3) 

and so 

f 2 R((3)- 1 > fiCiP)- 1 (7.4) 

We also note that R(j3) should be a better bound, for large f3 or small temperature, since 
it is linear in (3, whereas C(f$) is exponential in (3 and that this £(/3) could be understood 
better to explore the comparison of these two bounds. It also ionly depends on the potential 
<I> through the stability parameter B. 




Appendix A. Lagrange Inversion 

To find the inverse of the function se~ s as a power series, we can use the one-dimensional 
Lagrange inversion form from |Mer06] . We want to find a power series of s in terms of 
y = se~ s , which amounts to finding the inverse. We write the equation in the convenient 
form: 

a = ye s (A.l) 

For a formal power series s(y), implicitly defined by s — y<p(s), the Lagrange Inversion 
Formula can be stated as: 

[y"] S = V" 1 ]^)" (A.2) 
n 

where the notation [y n )f{y) is used to denote the coefficient of y n in the expression f(y). 
The function <fi on which we wish to perform Lagrange Inversion is <j>(s) — e s and the 
formula for the nth coefficient in the desired power series is: 

[y n ]s = - W 1 - 1 }^ 13 
n 

n [n — 1): n\ 

The power series we obtain is s = n n < V n as required. 

n>l 
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